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J> ■ Abstract 

qq \ We give an alternative description of the physical content of general relativity that does not 

require a Lorentz invariant spacetime. Instead, we find that gravity admits a dual description 
in terms of a theory where local size is irrelevant. The dual theory is invariant under foliation 
preserving 3-diffcomorphisms and 3D conformal transformations that preserve the 3-volumc 
(for the spatially compact case) . Locally, this symmetry is identical to that of Hofava-Lifshitz 
gravity in the high energy limit but our theory is equivalent to Einstein gravity. Specifically, 
we find that the solutions of general relativity, in a gauge where the spatial hypcrsurfaces have 
constant mean extrinsic curvature, can be mapped to solutions of a particular gauge fixing of the 
dual theory. Moreover, this duality is not accidental. We provide a general geometric picture 
for our procedure that allows us to trade foliation invariance for conformal invariance. The dual 
& ' theory provides a new proposal for the theory space of quantum gravity. 
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1 Introduction 

Conformal methods have a long and rich history in the study of gravity. Only two years after Ein- 
stein's seminal paper [1] on general relativity, Weyl had already proposed a theory [2, 3], invariant 
under conformal transformations of spacetime, as a way to unify gravity with electromagnetism. 
Unfortunately, no physically sensible theory with this symmetry has been proposed to date. At- 
tempts to implement 3D conformal invariance have met greater success. As early as 1959, Dirac [4] 
considered gauge fixing the spacetime foliations of general relativity using constraints that generate 
3D conformal transformations. These efforts culminated in York's discovery [5, 6] that the initial 
value problem for general relativity could be solved using conformal methods. Recently, spatial con- 
formal symmetry has played an important role in quantum gravity both through Hofava's proposal 
[7] for a UV completion of general relativity and Maldecena's AdS/CFT conjecture [8]. 

In this paper we show that, in gravity, there is a deep relationship between foliation invariance 
and 3D conformal invariance: they are different ways of expressing the same physical dynamics. 
The solutions of a particular gauge fixing of general relativity are equivalent to the solutions of a 
particular gauge fixing of a theory invariant under 3D conformal transformations. In the spatially 
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compact case, which we treat exclusively in this paper, the 3D conformal transformations must 
preserve the volume of space. This means that there exists an alternative way of understanding 
gravity that does not require a Lorentz invariant spacetime but relies instead on the relativity of 
local scale. We will demonstrate that the foliation invariance of general relativity can be traded 
for 3D local scale invariance and that the mechanism behind the trading has a simple geometric 
interpretation. Our result is inspired by York's solution to the initial value problem of general 
relativity. 

York began his 1973 paper [6] on the conformal approach to the initial value problem by 
stating: "An increasing amount of evidence shows that the true dynamical degrees of freedom of 
the gravitational field can be identified directly with the conformally invariant geometry of three- 
dimensional spacelike hypersurfaces embedded in spacetime." He continues: "the configuration 
space that emerges is not superspace (the space of Riemannian three-geometries) but 'conformal 
superspace'[the space of which each point is a conformal equivalence class of Riemannian three- 
geometries] x [the real line](ie, the time, T)." (Original parentheses.) Despite these bold claims, he 
does not show that the dynamics actually projects down to curves on conformal superspace nor 
does he provide an explanation for how the 4D diffeomorphism invariance of general relativity could 
be related to 3D conformal invariance. In this paper, we fill in these gaps. 

In a recent paper, Barbour and Murchadha [9] make use of an older result [10] to show 
that the dynamical solutions of general relativity on compact manifolds are uniquely determined 
by initial data on conformal superspace. This is possible provided the spacetime is foliable by 
spacelike hypersurfaces with constant mean extrinsic curvature (CMC). 1 They call this theory 
shape dynamics because it involves the evolution of local, scale-independent shapes. The first 
goal of our paper is to provide a rigorous geometric understanding of this result. We then extend 
this result by uncovering the precise relationship between general relativity and its dual theory - 
presumably Barbour and Murchadha's shape dynamics. 

We establish our result by introducing a Stiickelberg field 2 that is subsequently eliminated from 
the theory after performing a particular gauge fixing. This process and, in particular, the gauge 
fixing is inspired by the "best matching" used in [10]. For a detailed account of this procedure, 
that was first developed in [12], see [13] (or [14] for a canonical treatment). 

Our results both clarify and justify York's intuitive remarks regarding the configuration space of 
gravity: conformal superspace is not the reduced configuration space of general relativity but that 
of its dual theory. This fact could have important implications for quantum gravity. In the dual 
theory that we will develop, there is one global (non-local) Hamiltonian that generates the dynamics 
of the system on conformal superspace. The dual theory then suggests a natural candidate for the 
physical Hilbert space of quantum gravity: the set of states invariant under foliation preserving 
diffeomorphisms and volume preserving 3 conformal transformations (this also implies a similarly 
defined theory space). Given the many difficulties encountered in defining the physical Hilbert 
space of quantum general relativity, it may be useful to define quantum gravity as the quantization 
of the dynamics of shape. 

lr The restriction to CMC foliable spacetimes is weak since it includes the vast majority of physically interesting 
solutions to Einstein's equations while excluding many physically uninteresting solutions such as those that contain 
closed timelike curves or those that have exotic topologies. 

2 We use the term Stiickelberg field, because the mechanism is precisely the Stiickelberg's mechanism of restoring 
U(l)-gauge symmetry in massive electrodynamics [11]. 

3 They need only be volume preserving in the spatially compact case where the volume is well defined. 



3 



1.1 Outline 



The outline of the paper is as follows: we first illustrate our procedure for trading symmetries by 
applying it to a simple toy model. We then develop the general geometric picture that justifies our 
algorithm. Finally, we apply this algorithm to geometro dynamics and establish a dictionary that 
maps solutions of the dual theory to solutions of general relativity. In an appendix, we provide 
further mathematical details supporting the result on gravity. 

The toy model is a finite dimensional particle model constructed to exhibit many of the features 
of geometro dynamics. We've included it both for pedagogical purposes and because it provides an 
explicit example where all the equations can be worked out and the symmetry trading can be 
verified directly. However, this section can be skipped on a first reading if the reader is interested 
only in the main result. 



2 Duality in Toy Particle Model 

The toy model consists of n independent harmonic oscillators in d dimensions. The action is 4 

(1) 

The masses have been set to 1, for convenience. We will need to consider negative spring constants 
k < to avoid imaginary solutions. Unstable solutions will not worry us here since we are interested 
mainly in the symmetries of the toy model, which we do not consider a physical theory. 

We write the above action in a reparametrization invariant form by introducing the auxiliary 
fields Nj and the arbitrary label A 

where we have restricted to strictly negative k's. We choose to parametrize the trajectory of each 
particle independently and allow the time variables to be varied dynamically. The iVj's represent the 
A derivatives of the time variables and mimic the local lapse functions of geometrodynamics, whose 
terminology we will borrow. The total energy of each particle £i essentially specifies how much time 
should elapse for each particle in the interval Ai — Ao so that this theory is equivalent to n harmonic 
oscillators with different "local" time variables. We have introduced the local reparametrization 
invariance as a way of studying this feature of geometrodynamics in a finite dimensional toy model. 
It is not meant to have any direct physical significance. 

After performing the Legendre transform S = J ^A^™ =1 (n^qf — Ho(q,ir)^, where tt^ = to, 
the local reparametrization invariance leads to n first class constraints, xi- The total Hamiltonian 
is 

n 

H = Y, Nl Xl, where (3) 

7=1 

Xl = \tf-\k\q 2 i)-£l (4) 

4 Lowercase indices are spatial ranging from 1 to d and repeated indices are summed over. Uppercase indices label 
particle number and are only summed when explicitly shown. Also, v 2 = v a v b S a b for vectors and u 2 = u a UbS ab for 
covectors. 




(2) 
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and the phase space F(q, n) is equipped with the Poisson bracket 

{q%7r J b }=5 a b 5j. (5) 

Our goal, as explained below, will be to exchange all but one of the first class x/' s f° r new first 
class conformal constraints. 

2.1 Identify the Symmetries to Exchange 

Our stated objective is to exchange the local reparametrization invariance for conformal invariance. 
A naive choice would be to exchange all the %/'s f° r the symmetry, S, defined by identifying the 
configurations 

qj -> e+'(fi (6) 

as members of an equivalence class for some arbitrary scale parameters (pj. This choice, however, will 
lead to frozen dynamics because one global x must be left over to generate dynamical trajectories 
on phase space 5 . Thus, we need to leave at least one linear combination of x's unchanged. Because 
of this, our first-and most subjective-step is to identify a subset of S to exchange. 

A natural choice is the symmetry S/V, where V represents all configurations that share the 
same total moment of inertia. Since the total moment of inertia sets the global scale of the system, 
S/V identifies all locally rescaled g's that share the same global scale. This symmetry requires 
invariance under the transformation 

q a i -> e^q% (7) 

where <f> obeys the identity 

= 0. (8) 



The mean operator, is defined as (•) = — 2j?=i 'ii w here • signifies inclusion of the field we would 
like to take the mean of. We can impose the identity (8) explicitly by writing (f> in terms of the 
local scale factors (pi 

h = h- (<f>) • (9) 

Thus, we parametrize the group S/V redundantly by S. 
2.2 Introduce Stiickelberg Field 

The next step is to artificially introduce the symmetry S/V by applying the canonical transforma- 
tion t 

q a Uq a = Q a = q a 

iri -> t^i = Ui = e-Ki (10) 

generated by 

i 

We then enlarge the phase space T(q,ir) — > T e (q, ir, <fi, tt^) to include the canonically conjugate 
variables <pi and 7r^ that we use to construct <fi using (9). The non-zero Poisson bracket's introduced 
in T e are 

{fa,4} = 5j. (12) 



5 We draw the reader's attention to [15], which suggests that such a global x does not generate a gauge transfor- 
mation but a genuine physical evolution. 
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We can naturally extend t to a map t : T c — > T e by acting in the following way: t : (q, ir, (J), n^) i— > 
(t^qjt^ir, (j),^^). This gives the natural extension of functions, / £ C°°(r), on T to functions, 
f = fopot£ C°°(r c ), on T e induced by the canonical transformation (11): 

f(q, vr, cf), vr^) = fit^q, t^ir). 

Here, p : T e := T x T*5 — >• V is a canonical projection. From now on, we will abbreviate the 
transformation 6 above (p o t)* := T. 

The Poisson structure of T e is such that 

{Tf(q, tt), 4} = {Tf(q, tt), tt^? - (vr • g)} , (13) 

In the above, we have used the fact that n^qf — (ir ■ q) is invariant under T, as can be checked 
explicitly. Consistency between the dynamics and this Poisson structure requires the constraints 7 

C 1 = ttJ - (Tr^f - (vr • q)) « 0. (14) 

Given (13), the C s are trivially first class with respect to the transformed constraints T\i 

Xl^T Xl = l -(e-*Kl-e^\k\qi) (15) 

as they commute with any f(Q,U). This should be expected: the canonical transformation (11) 
did not change the theory. Since they act trivially on the functions extended into r e , the original 
theory is unmodified. But, by extending the phase space, we have added two phase space degrees 
of freedom per particle. The role of the C/'s is precisely to remove these degrees of freedom. The 
<j) field is essentially a Stueckelberg field for the conformal symmetry. This can only happen if the 
C/'s are first class. Consequently, the infinitesimal gauge transformations generated by the C/'s on 
the extended phase space variables 

SejQi = -5 IJ 0iqi 5 9j tii = 6 ij § i tt i 

5 6j <t>i = 5 IJ 8j b B ^\ = 0, (16) 

where 5oj- = {-,9 J Cj} (no summation), represents at this stage a completely artificial symmetry. 

Before writing down the remaining constraints, we note that one of the C/s is singled out by 
our construction. Averaging all the C/'s gives 

<^> ~ 0. (17) 

However, the action of this constraint is trivial for arbitrary functions / on the extended phase 
space r c . To see this, note that the gauge transformations 5q- = {■,6(ir ( p)} generated by (17) 
are trivial for the phase space variables q, tt, <fi, and tt^. The only transformation that is not 

automatically zero is 5g(f)j = 6 |<^/, (tt^) |, which can easily be seen to vanish. Since the </>'s only 

enter the theory through the <^>'s, we see that the constraints (17) act trivially on all quantities in 

6 This transformation is analogous to the restoration of [/(l)-gauge symmetry in massive electrodynamics through 
the indroduction of a Stiickelberg field (f>: The original photon field A M is replaced by a (7(l)-transformed + d^cj), 
where <p is not viewed as a gauge parameter but as a new field with its own gauge transformation property, such that 

2 

the transformed photon mass term -^-(A^ + d^^i^-fj, + d^cp) is gauge- invariant. 

7 These constraints could have been derived as primary constraints by introducing cf) before the Legendre transform 
through "best matching". See [13] or [14] for details. 
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our theory. This is a reflection of the fact that there is a redundancy in the <f> variables. Eq (17) is 
then an identity and not an additional constraint. This reflects the redundancy in the original (ft 
variables expressed in terms of 0's. 

It will be convenient to separate this identity from the remaining constraints by introducing 
the numbers af, where i = 1 . . . (n — 1), such that the set {(vr^) , Cj}, where Cj = ^ 7 a\ Cj, forms 
a linearly independent basis for the C/s. 8 We are then safe to drop the constraint (tt^) = from 
the theory. The remaining constraints are now 

T X i « 0, and d « 0. (18) 

where I = 1, . . . , n and i = 1, . . . , n — 1. 



2.3 Impose and Propagate Best— Matching Constraints 

In order to exchange symmetries, we impose the best-matching (also known as Mach constraints) 
constraints 

4 « 0. (19) 

These constraints are motivated by the desire make the theory background independent with respect 
to the C/V symmetry. For more details on how these constraints achieve this see, [16, 14]. One 
of the best-matching constraints, (71^) ~ 0, is first class but trivially satisfied. It can be dropped. 
The remaining 7r^'s can be represented in a particular basis as 71^ = djir^ and are second class with 
respect to the T^'s. To see this, consider the smearing functions f 1 . Then, the Poisson bracket 



n n 

i,j=i j=i 



f'TxiM**} = E [( a jf J W + \ k \ £)) - a jf J + \ k \ &] ( 2 °) 



is non-zero for general / (in the above, the mean is taken with respect to the j components). 
Nevertheless, the constraints 7r^ can be propagated by the Hamiltonian by using the preferred 
smearing f 1 = Nq that solves the n — 1 equations 

n 

£ cJjNi {\k\ + nj) = (a^N J (\k\ + . (21) 
J=l 

These lapse fixing equations lead to consistent equations of motion. 

We can now understand the reason for imposing the identity ^\ = 0. Were (f> unrestricted, 

the n x n matrix |Tx/,vr^| would have a trivial kernel A^q = 0. This uninteresting solution leads 

to frozen dynamics. Alternatively, the (n-l)xn matrix |Tx/,7r^| has a one dimensional kernel. 

Thus, (21) has a one parameter family of solutions parameterized by the global lapse. This leaves 
us, in the end, with a non-trivial time evolution for the system. 

The fact that we have a remaining global lapse indicates that we have a remaining linear 
combination of Tx/'s that is first class. This is the linear combination given by the preferred 
smearing A^q. We can split the T^/'s into a single first class constraint 

T Xi .c. = ^2 N o T Xi (22) 



8 Note that the ability to explicitly construct this basis is a feature of the finite dimensional model. This possibility 
sets the toy model apart from the situation in field theory, which is more subtle. 
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and the second class constraints 

T Xl = Ysti T Xi> ( 23 ) 
I 

where the (if are numbers chosen so that the above set of constraints is linearly independent. Note 
that, in the finite dimensional case, it is easy to ensure that this new basis is equivalent to the 
original set of Xi"' s - I n geometrodynamics, the situation is more subtle because we are dealing with 
continuous degrees of freedom. Using this splitting, we now have the first class constraints 

T X f.c. « 0, and Q « (24) 

as well as the second class constraints 

T Xi « 0, and ttJ « 0. (25) 

2.3.1 Excursion: Best Matching 

Since best-matching is just an inspiration for the construction in our paper, we restrict our expla- 
nation of best-matching to the minimum. 

In the construction of an action principle one is often faced with the problem that one can not 
directly parametrize the physical configuration space of a system. However, one can often describe 
the configurations space as the quotient of a redundant parametrizable space by the action of a 
gauge group, that one can parametrize as well. A prime example for this is geometrodynamics, 
where the distinct configurations, i.e. points in superspace, are not easily accessible, but can be 
described as diffeomorphism orbits in the space Riem, which consists of all 3-metrics. Each of these 
diffeomorphism orbits is a subspace of Riem. If we want to construct a geometrodynamical action, 
we have to construct it as an action for these subspaces. 

Let us now for simplicity assume that the action principle is given by the geodesic length 
of a curve in the redundant configuration space, for some given metric, and that the physical 
configurations are represented by non-intersecting subspaces of the redundant configuration space. 
The action principle asks for the minimal distance between two subspaces, so we look for the 
shortest (best-matching) geodesic that starts at one and ends at the other. The action thus 
gives two conditions (1) the curve has to be a geodesic and (2) it has to intersect both subspaces 
orthogonally. These two conditions generalize to (1) the curve has to extremize the action and (2) 
the conjugate momenta have to be orthogonal 9 . One can straightforwardly see that if 7 parametrizes 
the gauge group (and hence the orbit manifold) that this condition on the conjugate momenta is 
equivalent to demanding that the best-matching condition 

dL 

holds. One can furthermore straightforwardly see that for the 3-diffeomorphism group the condition 
one arrives at is the momentum constraint. For more details on this derivation, see [17]. 

2.4 Eliminate Second Class Constraints 

To see that this procedure has succeeded in exchanging symmetries, we can eliminate the second 
class constraints by defining the Dirac bracket 

n-l 

{•> -lob = {■, "}pb + £ {-'41 C l i T Xi, -} - {;T X i} C) {ttJ, ••} , (26) 



9 As a one-form, it in fact has to have the submanifold as its milliner. 
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where 

C5 = {4,T Xj -}- 1 . (27) 

The existence and uniqueness of the Dirac bracket depend on the invertibility of C, which, in turn, 
depends on the existence and uniqueness of the lapse fixing equations (21). In this toy model, 
(21) are just arbitrary linear algebraic equations. Thus, C will generally have a non-vanishing 
determinant. 

The advantage of using the Dirac bracket is that the second class constraints tt % ^ and T\i become 
first class and can be applied strongly and eliminated from the theory. This procedure is greatly 
simplified because one of our second class constraints, namely tt 1 ^ = 0, is proportional to a phase 
space variable. Thus, the method discussed by Dirac in [4] can be applied to eliminate the second 
class constraints. For 7r^, this involves setting tt\ = everywhere in the action. For the T%i's, this 
involves treating T\i as an equation for <fi. The solution, <pQ(q,n), of 

Txil^i = (28) 

is then inserted back into the Hamiltonian. This leaves us with the dual Hamiltonian 

n— 1 n 

#duai = NT X fA<h) + E A * E d i (9>a " (Q ■ T» , (29) 

i=l 1=1 

where N and A* are Lagrange multipliers. We can write this in a more convenient form by using 
the redundant constraints Dj = q^TT^ — (q ■ tt) 

n 

tfdual = NTxLcXM + E A ' D I ( 30 ) 

1=1 

remembering that one of the Dj's is trivially satisfied. 

The remaining Hamiltonian is now dependent only on functions of the original phase space 
r. Furthermore, the Dirac bracket between the dual Hamiltonian and any functions /(<?,7r) on T 
reduces to the standard Poisson bracket. This can be seen by noting that the extra terms in the 
Dirac bracket {f(q, tt), H^ USi i} Dh only contain a piece proportional to {f(q, tt), tt^}, which is zero, 

and a piece proportional to, {ttJ,, #dual}> which is weakly zero because -ffdual is first class. Thus, 
this standard procedure eliminates <f> and tt^ from the theory. 

The theory defined by the dual Hamiltonian (30) has the required symmetries: it is invariant 
under global reparametrizations generated by Xf.c. |^ an d it is invariant under scale transformations 
that preserve the moment of inertia of the system. This last invariance can be seen by noting that 
the first class constraints Dj generate exactly the symmetry (7). 

Before leaving this section we will give two explicit examples of how one can construct Txt.cX^o) 
explicitly. In the first example, we solve the second class T^j's for 0q exactly and then plug the 
solution into the definition of Txt.c- This immediately gives the general solution for n particles. 
Unfortunately, this procedure cannot be implemented explicitly in geometrodynamics because the 
analogue of Txi((fto) = is a differential equation for <f>Q. There is, however, a slightly simpler 
procedure that one can follow for constructing Txf. c .(0o) provided one can find initial data that 
satisfy the initial value constraints Xi of the original theory. If such initial data can be found, then 
(f>Q = is a solution for these initial conditions. Because the theory is consistent, the initial value 
constraints will be propagated by the equations of motion. Thus, 0o = for all time. We can then 
use this special solution and compute Xf.c. using its definition. This alternative method is given in 
the second example for the n = 3 case and agrees with the more general approach. 
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2.4.1 General n 



For simplicity, we take £j = 0. This will lead to transparent equations, which we value more in the 
toy model than physical relevance. In the particle model, it is possible to find <pQ explicitly for an 
arbitrary number of particles, n, by solving Txi{4>o) = 0. It is easiest to solve for <pQ directly. The 
solution is 10 



(31) 

We can choose a simple basis for the second class x/ s: {.XiV 7^ i}> f° r some arbitrary j. The first 
class Hamiltonian is then 

Xt.c=Nlxj(4&). (32) 
Using the identity <jj = — V i we can rewrite this as 

Xf.c. oc exp ^2 ^ $ j tt] + exp ^-2 £ % \ \k\ q). (33) 
Inserting (31), the first class x becomes 

n n 

i=i i=i 

It can be readily checked that (34) is invariant under the symmetry 

qi — > e^'qi 717 — > e _( "7i7, (35) 

where, (j) 1 = (f) 1 — {(f)). 

In geometrodynamics, solving for eft explicitly will not be possible because it will require the 
inversion of a partial differential equation. However, if one is only interested in mapping solutions 
from one side of the duality to the other, then it is sufficient to find initial data that satisfy the 
X/'s with (f> = and use the preferred lapse to construct the global Hamiltonian. This procedure is 
outlined in the next example. 

2.4.2 Example: n = 3 

We analyse the 3 particle case because this involves subtleties that do not arise in the 2 particle 
case but must be dealt with in the general n particle case. 

If we choose the basis {vr^} = { 7^,71"^}, then the lapse fixing equations (21) with = have 
the unique solution 

N ° ~ N ° [\k\ q\ + N °- N ° [\k\ q\ + A ) ■ (36) 

Inserting this solution into Xf. c .|^ = o = J2i NqXi, we nn d 

Xf.c oc (ttI^I - |fc| 3 qlqlqf) - 2xiX2Xz- (37) 



Negative spring constants k < are required to give real solutions to this equation. 
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The consequence of setting (ft 1 = instead of using it to solve Txi{4>q) = is that we pick up a 
term that is proportional to xi an d X2, which are strongly zero. For this reason, we need to ensure 
that we have initial data that solve the initial value constraints (xi = and xi = in this case) 
for eft = in order to arrive at the correct global Hamiltonian. Assuming that such data have been 
found, the last term is zero and Xf.c. reduces to 



/ 2 2 2 I ; 1 3 2 2 2 \ 

Xf.c oc (jr^TTg - \k\ <7ig 2 <?3 ) ■ 
This agrees with our general result from last section. 

2.5 Construct Explicit Dictionary 

We can now compare our starting Hamiltonian with that of the dual theory: 

i 

#dual = N Xf.c. I to + E A " W«? " fa • «)) 



(38) 



(39) 



N 



i=i 



i=i 



+ ^A J {nlq? - (7r 7 • qj)) . 



(40) 



The first theory is locally reparametrization invariant while the second has local scale invariance 
and is only globally reparametrization invariant. Note the highly non-local nature of the global 
Hamiltonian of the dual theory. 

The equations of motion of the original theory are 



qf = N\l8 ah 



it I a = N I \k\q b I 5 ab . 



Those of the dual theory are 



Qi 



2N U 



7Ti 



„I tab , n ~ 
7T b d H 



n 



k iq a j 



2N]J \k\q 2 j 



I A: | q%b + ^^A/vrf . 
n 



(41) 



(42) 



(43) 



We can now read off the choice of Lagrange multipliers for which the equations of motion are 
equivalent 



N I = 2NH7r 2 j^2NH\k\q 2 j 
A 1 = 0. 

In the second equality of (44), the xi' s have been used. 



(44) 
(45) 



3 General Symmetry Trading Algorithm 

Based on the above description of the procedure, we can sketch out a general algorithm for trading 
the first class symmetry x f° r D: 
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1. 



Artificially introduce the symmetry to be gained in the exchange by performing a canonical 
transformation of the schematic form: 



F((j>)=qexp((i>)Il. (46) 

This extends the phase space of the theory to include the Stueckelberg field <j> and its conjugate 
momentum tt^ and introduces constraints of the form 

C = D - vr^. (47) 



2. Impose the best-matching conditions 

^ « 0. (48) 

If they are first class, then the procedure just "gauges" a global symmetry. If they are second 
class with respect to the transformed first class constraints T\ of the original theory, then 
they can be exploited for the exchange if the Lagrange multipliers, iV , can be fixed uniquely 
such that 

{T X (N°),7T <p }=0. (49) 

When this is possible, the tt^s can be treated as special gauge fixing conditions for x- It is 
possible, as in our case, that there will be a part of T\ that is still first class with respect to 
7r^. This should be separated from the purely second class part. 

3. Define the Dirac bracket to eliminate the second class constraints. This can be done provided 
the operator {Tx,Tr&} has an inverse (which also means that N° exists and is unique). The 
second class constraints can be solved by setting tt^ = everywhere in the Hamiltonian and 
by setting cj> = </>o such that Tx(4>o) = 0- The second class condition and the implicit function 
theorem guarantee that this can be done. When this is done, the Dirac bracket reduces to 
the Poisson bracket on the remaining phase space variables. The 7t^, terms drop out of the 
C's and the x' s have been traded for .D's. 

4. To check consistency, construct the explicit dictionary by reading off the gauges in both 
theories that lead to equivalent equations of motion. 

Remark: We first pass from an original gauge theory on phase space T to a theory on a larger 
phase space T e , containing the Stiickelberg field 4> as a new physical field and not a gauge parameter. 
The new theory is however precisely the original if one imposes the C constraints and chooses a 
gauge such that = 0. 

Guided by this basic algorithm, we will now construct a formal geometric picture to illustrate 
why this algorithm works. 



3.1 Geometric Picture 

Let us now examine the mathematical structure behind the trading of gauge symmetries that we 
encountered by constructing the dual theory. We start out with a gauge theory (r, {.,.}, H, {xi}iez), 
where T is the phase space of the theory supporting the Poisson-bracket {., .}, the Hamiltonian H, 
and the constraints Xi- We demand that the constraints are first class, i.e. for all i,j £ I there 
exist phase space functions as well as s.t. 

{XiiXj} = Hk&xfijXk ^q-j 

{H,Xi} = Y.k&i u iXk- 
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The constraints define a subspace 11 C = {x S T : Xi(^) = Vi G X} and the first line of equation (50) 
implies that the derivations {xi, •} are tangent to C, while the second line implies that the Hamilton 
vector field {H, .} is tangent to C. The {xi, •} generate a group G of gauge-transformations on C 
that lets us identify C through an isomorphism i : E — )• C with a bundle 12 E over C/G. The fibres 
of E are gauge orbits of a point x £ C and thus isomorphic to G/Iso(x). According to Dirac, we 
identify the fibres with physical states and observe that the total Hamiltonian Htot — H + Yliei ^Xi 
depends on a set of undetermined Lagrange multipliers {Aj}j e j. Fixing a gauge means to find a 
section a in E. 

This is most easily done by imposing a set of gauge-fixing conditions 13 {(/>i}j g j such that the 
intersection of Q = {x G T : 4>i(x) = Vi € 1} with C coincides with i(cr). To preserve i(cr) under 
time evolution we have to solve {H tot , 4>i}\ c = \/i £ I for the Lagrange-multipliers A 1 = Xi. The 
gauge-fixed Hamiltonian is 

//,/: " • »:,\, (51) 

and the equations of motion are generated by {H g f, .}, while the initial value problem is to find 
data on z(cr). 

Given a gauge theory (H, {xi}iei) on a phase space (r, {., .}) we define a dual gauge theory as a 
gauge theory (Hj, {pj}jej) on (r, {., .}) if and only if there exists a gauge fixing in the two theories 
such that the initial value problem and the equations of motion of both theories are identical. 

One particular way to construct a dual theory is to use gauge-fixing conditions {4>i}i£i defining 
Q such that for all i,j G I there exist phase space functions gfj satisfying the integrability condition 

{&,0j} = X)4&- (52) 

fcgZ 

The {4>i, .} thus generate a group EI of transformations on Q that lets us identify Q with a bundle 
F over Q/M using an isomorphism j : F — >■ Q. The fibres of F at x £ Q are M/Iso(x). The 
Hamiltonian of the dual theory has to satisfy 

{Hd,-}\ a = {Hgf,.}\ a , . 

{H d ,^}\ g = 0. ^ 

These conditions as well as the integrability condition can be fulfilled by construction using the 
canonical Stuckelberg formalism to implement symmetry under an Abelian group H with the sub- 
sequent elimination of the Stuckelberg field by substituting it with the solution to the constraint 
equations in the original theory. The nontrivial conditions for the Stuckelberg formalism to yield 
the anticipated trading of gauge symmetries are: 1) that the generators of the group furnish a 
gauge-fixing for the original gauge symmetry and 2) that the transformed constraint equations 
admit a solution in terms of the Stuckelberg field. In the particle model, the first condition is 
satisfied by the invertibility of (21) in terms of Nq while the second condition is guaranteed by the 
invertibility of Txi{4>o) = in terms of <f>o. 



Although we assert "space" we abstain from topologizing T or any of its subsets in this paper. 
12 This bundle is, in general, not a fibre bundle, since different points x in C generally have different isotropy groups 
Iso(x). 

13t 



In general the index set for (f>i could be different from I. 
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4 Gravity as a Conformally Invariant Theory 



We now apply the general procedure to geometrodynamics. Before applying the procedure, we 
establish notation and review canonical general relativity. We treat in detail the spatially compact 
case leaving the asymptotically flat case for future investigations. For more details on the mathe- 
matical structure of the dualization and a more detailed proof of the dictionary, see Appendix (A). 

We start with the ADM formulation of general relativity on a compact spatial manifold E 
without boundary (and if confusion arises we will assume the topology of E to be 5 3 ). The 
phase space T is coordinized by 3-metrics g, represented locally by a symmetric 2-tensor g : 
x ' — y g ab (x)dx a dx b , and its conjugate momentum density 7r, represented locally by a symmetric 
2-cotensor tt : x \— > 7r ab (x)d a d b of density weight 1. Given a symmetric 2-cotensor density, F, and 
a symmetric 2-tensor / we denote the smearing by 

?ab/„\„ — a „( f\ / ^3„„a6 



F(g) := / d 6 xF a \x)g ab {x) and vr(/) := / d 6 xir ab (s) f ab (x) . (54) 

We will not explicitly state differentiability conditions for (g, tt) or details about the Banach space 
we use to model T, we just assume existence of suitable structures to sustain our construction. The 
nonvanishing Poisson bracket is 

{F(g),7r(f)} = F(f) := / d 3 xF ab (x)f ab (x), (55) 

and the Hamiltonian is 

H(N,0= [ d 3 x(N(x)S(x)+C(x)H a (x)), (56) 

JT, 

where the Lagrange mutipliers TV and £ a denote the lapse and shift respectively. The constraints 
are 

S{x) = K ab {x)^fM^ d (x) - y/\${x)R\g\{x) 

VlslW (57) 
H a (x) = -2g ac (x)D b Tr bc (x), 

where D denotes the covariant derivative w.r.t. g, G abc d denotes the inverse supermetric and R[g] 
the curvature scalar. Denoting smearings as C(f) = J E d 3 xC(x)f(x) and C[v) := J T .d 3 xv a (x)C a (x), 
we obtain Dirac's hypersurface-deformation algebra 

{H(u),H(v)} = H([u,v\) 

{H(v),S(f)} = S(v(f)) (58) 
{5(A), S(f 2 )} = #(iV(A,A)), 

where [.,.] denotes the Lie-bracket of vector fields, so the first line simply states that the H(x) 
furnish a representation of the Lie-algebra of vector fields on E and where -/V a (AjA) : % l— >• 
<f\x)(h{x)hAx)-hAx)f 2 {x)). 



4.1 Dualization 

4.1.1 Identify Symmetries to Exchange 

Let us spell out the symmetry to be gained in exchange for foliation invariance. Naively, we would 
trade S{x) for constraints that generate general conformal transformations of g. However, as in 
the toy model, trading all such symmetries would lead to frozen dynamics. One global constraint 
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must be left over to generate global reparametrization invariance. In analogy to the toy model, we 
restrict to conformal transformations that do not change the global scale. In geometro dynamics, the 
analogue of the moment of inertia is the total 3-volume. Thus, the desired symmetry is explicitly 
constructed in the following way: 

Let C denote the group of conformal transformations on S and parametrize its elements by 
scalars 4> : £ — > R acting as 

' 9ab(x) -> e 4 ^g ab (x) 



ir ab (x) 



-4<f>(x) ^ab / x y 



(59) 



The factor of 4 is conventional in 3D. 14 It is chosen so that the scalar curvature R[g] transforms in a 
simple way under the above tranformation. Consider the one-parameter subgroup V parametrized 
by homogeneous <f> : x — > a. Notice that V is normal, because C is Abelian, so we can construct the 
quotient C/V by building equivalence classes w.r.t. the relation 



C 3 (f> ~ <f>' iff 3a G V, s.t. 



+ a. 



(60) 



Given a metric g on £ and <j) G C we can find the unique representative in the equivalence class 
[4>]~ € C/V that leaves V g = J s d 3 x^/\g\(x) invariant using the map 



- ln(e 



6<M 



///- 



(61) 



where we define the mean (f) g := y J^d?xy f \g\{x)f{x) for a scalar / : S 



The map 



allows us to parametrize C/V by scalars (p. Note that (j) g can be written more transparently by 
observing that it is chosen so that the volume element of the conformally transformed metric is 
equal to 

M 



1.91 



Thus, the conformal transformation 



(62) 



9ab{x) ->■ exp (40 (a;) ] g a b(x) 



leaves 



V„ 



d xy\g\ 



d 3 x ' 



(>4> 



\g\ = v 9 



(63) 



(64) 



invariant (we will often suppress the subscirpt g in <p g for convenience). C/V is then precisely the 
symmetry we want to obtain in exchange for foliation invariance. 



4.1.2 Introduce Stiickelberg Field 

The next step is to perform a canonical transformation that artificially introduces the symmetry 
C/V by inserting a Stiickelberg field (f). If we define the generating function 

F[4>) := J d 3 xg ab (x) exp (40(s)) U ab (x), (65) 
14 In dimension d, the conventional factor is -rro- 
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then the canonical transformation t generated by F[(f>] is: 



9ab(x) 
TT ab (x) 



tg a b(x) 
tir ab (x) 



exp [4<t>(x)J g a b(x) 

(V 6 (z) - i(vr). 



1 



9\{x, 



(66) 



This is indeed a volume preserving conformal transformation. 

We now enlarge the phase space T — > T e = V x T*(C) with the canonical 'product Poisson 
bracket' and parametrize conformal transformations by scalar functions <p and their conjugate 
momentum densities by scalar densities tt^. As in the particle case, to extend functions from T to 
r e we define T as (po t)*, where p denotes the canonical projection p :T x T*(C) = T e — > T. Note 
that p* is also a canonical transformation in the sense that the Poisson bracket of its image in T e 
coincides with the Poisson brackets in T. 

The scalar constraints transform in the following way under the map T: 



S(x) -> TS(x) 



-6<t> 



7T Tab 



IT 

2 



M (1 

6 V 



=60 ^ 



R[g}-8 [D 2 <f>+ (D(f>f 



(67) 



where v 2 = g a bV a v b for any vector v a . To avoid technical difficulties arising from the proper 
treatment of the diffeomorphism constraint, we relax it for the time being and verify after the 
dualization is completed that it can be consistently reintroduced into the theory. 
The quantity 

V = n(x)-V\g~\(x)(7r) g (68) 

will be the first class constraint left over in the dual theory. We pause for a moment to note its 
important properties, which can be verified by straightforward calculations. First, T> is invariant 
under the canonical transformation (66). Second, it generates infinitesimal volume preserving 
conformal transformations. This can be seen by noting that 



6 e g a b(x) = mx)-(49))g ab (x) 



5 e ^ b {x) = (-46(x) + (40)) ( ^ b {x) - - (ir) g g a \x)VW\(x) 



(69) 



where Sg = {-,X>(40)}, is the infinitesimal form of (66). This is the key property that we require of 
the dual theory. 

For the moment, the relevance of the above observation is that it allows us to show that the 
canonical transformation has introduced a new constraint. From (69), it is simple to verify that 
the action of T> on arbitrary smooth functions / of T transformed by T is equivalent to the action 

Of 7T^ 

{17, vr^(x)} = {T/, tt(x) - yfig\(z)(n) g }. (70) 



Thus, the quantity 



C(x) = 4 (n(x) - V^T(x)(7r) ff ) - ir^x) 



(71) 



acts trivially on the image of T. In order for this to be consistent with the dynamics, C(x) should 
be added to the Hamiltonian as a constraint using a Lagrange multiplier. This is analogous to the 
toy model: the canonical transformation has both introduced auxiliary quantities (f>(x) and n Ax) 
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and corresponding constraints C(x). Because of (70), C(x) is trivially first class with respect to 
any transformed quantities. Thus, 

{C(f 1 ),TS(f 2 )} = (72) 
{c(f),TH(v)}=0 (73) 

for the smearing functions /i,/2,/, and v. Note that C(x) could have been derived in a more 
conventional way as a primary constraint by introducing the Stiickelberg field before performing 
the Legendre transform. This equivalent approach was followed in [10]. 

We note one final property of C{x). Just as in the toy model, one of the C(x)'s is an identity 
and not an independent constraint. This can be seen by noting that 

/ d 3 xC(x) = [ d 3 xTr^x) (74) 
Jt, Jt, 

but the Poisson bracket of J d 3 xir ( j > with the variables gij, ir^ , <p, and tt^ is identically zero (the 

only non-trivial calculation is | /)((/>), J* s d?x 7r^(x)| = for arbitrary smearings p). This is a 

consequence of restricting to C/V instead of just C. It means that / afx C can be removed from 
the theory without affecting the theory in any way. The result of this is that we will be left with a 
global first class scalar constraint at the end of the procedure. 



4.1.3 Impose and Propagate Best Matching Constraints 



We now carry out a gauge fixing inspired by best-matching [10, 16, 14]. This involves imposing 
the constraints 

Tty = 0. (75) 

One of these constraints: J ci 3 x7r^, is trivial as noted above. However, singling out this constraint 
explicitly is more subtle now than it was in the finite dimensional toy model. For simplicity, 
we will work with this redundant parametrization of the constraints, keeping in mind that it is 
over-complete by one equation. 

Imposing 7r^(x) = turns all but one of the original scalar constraints into second class con- 
straints, as can be observed from the Poisson bracket 



{TS(N)^(\)} = / dx\(x) 



F N 



\g\ix e 



6kx) (F N ) g 



(76) 



where we smeared ir^ with a scalar A and Fjv is given by 



F N = 8g ab D a (e^D b N^j s/\g~\ - 8Ne 2 ^ 



R[g] - 8 [D 2 ^ + (D^f 

2Ne^^/\g\(*) 2 9 ~ N [6T5 + 2^> 3 (C + vr )] . (77) 



The last term in Fn is weakly zero. The consistency of the dynamics requires that we fix the 
Lagrange multipliers to satisfy the equation 



{TS(N ),^(x)} 



0. 



(78) 



The trivial solution Nq = yields "frozen dynamics." Fortunately, there is a unique solution 
Nq((/), gij, it 13 ) to the lapse fixing equation 



N 



-40 



R[g}-%{D 2 ^ + {D^ 2 )] +- A ^)\ 



e- & ha b D a (e 2 *D b N) = {£) 



(79) 
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that is equivalent to the first line of (78). In (79), C is defined as e^\/\g\ times the left-hand-side 
of (79). 

In terms of the transformed variables G a b = Tg^ and Ii ah = Tir ab , (79) takes the simple form 

(IB- n\ I ( r , (II) r o 



R[G] + — D J N = (^ yR[G] + - D z j N J . (80) 

This is the same lapse fixing equation obtained in [10] by a similar argument using the Lagrangian 
formalism. To understand why (80) has a one parameter family of solutions (and, thus, depends 
only on </>) note that any solution, Nq, to 

R[G] + ^-dAn = -c 2 , (81) 



for an arbitrary constant c, is also a solution to (80). Eq (81) is the well-known lapse fixing equation 
for CMC foliations. It is known to have unique positive solutions[18, 5, 19]. This fact is vital for 
our approach to work as it allows us to construct the Hamiltonian of the dual theory. 

The one parameter freedom in Nq is exactly what we expect from the redundancy of the 7T0(x)'s. 
We are one global equation short of completely fixing the N(xys. Thus, we should be left with a 
single first class linear combination of the T5(x)'s. This is precisely the gauge fixed Hamiltonian 



J d 3 xN (x)TS(x) = TS(N ), (82) 



where both Nq(x) and TS(x) are functionals of g,ir,(fr. H g { is constructed to be first class with 
respect to tt^. We can now state the important result that the extended system with the constraint 
= is consistent as merely a (partial) gauge fixing of the original system. 



4.1.4 Eliminate Second Class Constraints 

Let us collect and classify the constraints of our theory to recap what we have done. We started 
with the first class constraints TS(x) and then imposed the conditions n^(x) = 0. The former 
are second class with respect to the latter but, since J d 3 XTr ( f ) plays no role in the theory, there is 
a corresponding constraint H g { that is still first class. We can split the first class H g { from the 
remaining second class TS'(x) , s by defining the variable TS(x) = TS(x) — H g {. This leaves us with 
the first class constraints 

Hg w C(x) w (83) 

as well as the second class constraints 

TS(x) tt^x) w (84) 

in analogy with equations (24) and (25) of the toy model. Note that, for the moment, we are still 
relaxing the diffeomorphism constraints. In Sec (4.1.3), we showed that iVo is consistent gauge 
choice for general relativity. However, we should now construct the Dirac bracket to explore the 
full structure of the theory. 

First, notice that since {Hgf,^^} w we know that {TS(x),TT ( f ) (y)} w {TS(x),n ( j > (y)}. This 
leads to the weak equality 

J G(x, x'){TS(x'), ir ( / ) (y)}d 3 x' « J G(x, x'){TS(x'), ^(y)}d 3 x' = 5(x, y) 
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provided the Green's function, G, for the differential operator acting on N in (80) exists. This is 
guaranteed because the existence and uniqueness of solutions to the CMC lapse fixing equation, 
for suitable initial data and boundary conditions, implies existence of the Green's function for the 
respective initial data and boundary value problem. 15 The formal existence of this Green's function 
is all that we will need for the remainder of the paper. 

In terms of this Green's function, the Dirac bracket is defined as: 



{/i,/ 2 }d := {fij 2 }- J d*xd?y{h^{x)}G{x,y){TS{y),h} 

+ / d s xd 3 y {f u TS(x)}G(x, y){^(y), f 2 }, 



(85) 



for arbitrary functions f\ and f 2 - By construction, the Dirac bracket between tta or TS and any 
phase space function is strongly equal to zero. We can then eliminate these constraints from the 
theory by following Dirac's algorithm [4]. First, we set tt^ = everywhere that it appears; then 
we find (j)Q such that TS((j)o) = and insert </>o into all constraints. Note that it is (j), and not <fi, 
that we can explicitly solve for. This is analogous to what happens in the toy model and arises 
because of the redundancy in the tt^x). When this is done, the system will have been reduced to 
the original phase space T. The Dirac bracket between functions / of T and the remaining first 
class constraints is weakly equivalent to the Poisson bracket. This is true, just as in the toy model, 
because the extra terms in the Dirac bracket are either proportional to {/, ^(x)}, which is zero, 
or {^(y), H g {} and {^(y), n ( j ) (x)}, which are both weakly zero. Finally, because we are simply 
inputting the solution of a constraint back into the Hamiltonian, the equations of motion on the 
contraint surface will remain unchanged. Thus, T^ is still effectively a canonical transformation. 
For more details, see Appendix (A) and, in particular, the steps leading to Eq (110). 

The last step is to verify that TS((po) = can be solved for (pQ. We can simplify (67) using the 
strong equations tt^ = and by taking linear combinations of the first class constraints C(x). This 
leads to the equivalent constraint 



-60 



(vr) 



0, 



(86) 



where a ab = 7r ab - ± {ir) g g ab ^\g\ and R[g] = R[g] - 8(D 2 (j) + (D(f)) 2 ). Eq (86) is the Lichnerowicz- 
York (LY) equation used for solving the initial value problem of general relativity. Its existence 
and uniqueness properties have been extensively studied. 16 It is known to have unique solutions 
when a ab is transverse and traceless. Fortunately, these are exactly the conditions required by the 
diffeomorphism and conformal constraints respectively. The formal invertibility of this equation is 
the second vital requirement for our procedure. Without this, we would not be able to prove the 
existence of the dual theory. However, given that we can solve (86) for 4>q for specified boundary 
and initial data, we arrive at the dual Hamiltonian 



H' dual = MH gf [M + / d*x\(x) (tt(x) - <7T» 
E 



(87) 



where N is a spatially constant Lagrange multiplier representing the remaining global lapse of the 
theory. We can now reinsert the diffeomorphism constraint. This gives the final Hamiltonian 



F dual =AAtf ff/ k/> ]+ / d 3 x \(x)V{x)+e{x)T io H a (x 



15 This is implied by the existence of a Green function for the so called Lichnerowicz Laplacian D + R, an operator 
that can be put into the form of a Hodge Laplacian dS + Sd. 



id 



See [19] or, for the specific context given here, see [20]. 
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using the definition (68). Note that we must use the transformed H a {x) evaluated at cpo (recall 
that T> is invariant under T). As shown above, the usual Poisson bracket over T can be used to 
determine the evolution of the system. 

4.1.5 Construct Explicit Dictionary 

We can verify that the dual Hamiltonian has the required properties. First, we check that all 
the constraints are first class. The gauge fixed Hamiltonian was constructed to be first class with 
respect to the conformal constraints T>. This can be seen by observing that the T5'(x)'s are first 
class with respect to the C(a;)'s and that the C(x)'s are equal to the V(x) 7 s when 7r^(x) = 0. The 
HaS are easily seen to be first class with H g { because they are first class with respect to the original 
S(xYs and because T is a canonical transformation. The 7r^'s commute with themselves because 
they are ultralocal canonical variables. Lastly, one can directly verify that 

[T io H{v),VU)} = T^ Q [H(v),V(f)} = D{C v f) « 0, (89) 

where v and / are smear ings. 

Secondly, the dual theory is indeed invariant under volume preserving conformal transforma- 
tions. For this, recall that the T>'s generate the infinitesimal form of (66) according to (69). Fur- 
thermore, the theory is also invariant under 3D diffeomorphisms generated by the H^s. The theory 
is not, however, invariant under 4D diffeomorphisms. The diffeomorphism invariance is only within 
the spatial hypersurfaces and is, thus, foliation preserving. This means, in particular, that the 
theory is not Lorentz invariant because it is not invariant under boosts. 

Finally, there is a gauge in which the equations of motion of the two theories are equivalent. 
Compare the ADM Hamiltonian to that of the dual theory. 

H ADM = J d 3 x (N(x)S(x) + e(x)H a (x)) 

\{x)V{x)+e{x)T i: H a {x) 
= J d 3 x (NN (x, 4> )T^S(x) + \{x)V{x) + £ a (x)T^ a (x)) . (90) 

Because T is a canonical transformation, the equations of motion of both theories take the same 
form in the gauges 

N(x) =NN (x,fo) 

X(x) = 0. (91) 

To map the solutions from one side of the duality to the other, we simply need to use the explicit 
function <pQ and the map T^. This completes the dictionary. 

This dictionary is particularly straightforward to use if one can find initial data for the ADM 
Hamiltonian that satisfies the initial value constraints. In this case cpo = 0. Because both theories 
are first class, this condition will be propagated and the solutions of each theory are equal in the 
gauges (91). To find suitable initial data, we still must solve the LY equation. However, for the 
purpose of using the dictionary in the classical theory, we only need to solve this for a single point 
on phase space. 17 

Given the existence of the above dictionary, we arrive at the following proposition: 

17 To prove that the dual theory actually exists and to study the quantum theory, we still must be able to solve the 
LY equation over all of phase space. 



ffdual =NH gf [<j)o] + / d 3 X 
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Proposition 1 The theory with total Hamiltonian (88) is a gauge theory of foliation preserving 
3-diffeomorphisms and volume preserving 3D conformal transformations. In the gauge A = 0, this 
dynamical system has the same trajectories as general relativity in CMC gauge. 

5 Conclusions and Outlook 

Let us summarize the main result of this paper: we showed that there is a duality between gen- 
eral relativity and a gauge theory of 3D conformal diffeomorphisms that preserve the total spatial 
volume. Such a theory does not admit a local scale and is, therefore, a theory of the dynamics of 
the intrinsic "shape" of objects. We identify this theory with Barbour and Murchadha's shape 
dynamics. The classical duality means that there exists a distinguished gauge in both theories such 
that the gauged trajectories of one theory coincide precisely with the gauged trajectories of the 
other. In other words, one can trade the gauge symmetry of general relativity (spatial diffeomor- 
phisms and local refoliations) for the gauge symmetry of shape dynamics (spatial diffeomorphisms, 
volume preserving conformal transformations, and global refoliations) without changing the phys- 
ical content of the theory. The Hamiltonian of shape dynamics is nonlocal and difficult to write 
explicitly. However, explicit expresssions can be obtained, for example, by applying our procedure 
to the torus universe in 2+1 dimensions. One can draw a number of physical conclusions from this 
duality: 

1. Physical degrees of freedom of general relativity: The physical degrees of freedom of general 
relativity and shape dynamics are identical. This confirms earlier results [6, 10, 9] that the 
physical degrees of freedom of general relativity can be labeled by conformal 3-geometries 
and their momenta. 

2. Theory space of quantum gravity: The only input into the asymptotic safety approach to 
quantum gravity is the theory space. While the theory space of general relativity is usually 
identified with 4-diffeomorphism invariant metric theories, the theory space of shape dynamics 
is 3D conformal diffeomorphism invariant metric theories. This difference in theory space 
may lead to a different fixed point structure and different resolutions to the ghost problem of 
quantum gravity. Furthermore, the theory space of shape dynamics is locally identical to that 
of Hofava-Lifshitz gravity. This provides a possible link between Hof ava-Lifshitz gravity and 
shape dynamics. 

3. Loop Quantum Shape Dynamics: Following the Loop Quantization program, one can quantize 
shape dynamics. Observing that scale information in Loop Quantum Gravity is encoded 
in the spin quantum numbers on the edges and the intertwiner quantum numbers on the 
vertices of spin network functions, one may speculate that the physical Hilbert space of Loop 
Quantum Shape Dynamics is spanned by unlabeled graphs carrying a total volume quantum 
number. With this proposal, one should study Hamiltonians on this physical Hilbert space 
that propagate a CMC condition. 

4. New possibilities for quantum geometrodynamics: Since all constraints of shape dynam- 
ics generate transformations that preserve the configuration space of all 3-metrics, one can 
implement the gauge transformations as pullbacks on this configuration space. The corre- 
sponding operator algebra is naturally represented on a kinematic Hilbert space in which each 
distinct 3-metric defines a normalized geometry eigenstate orthogonal to all other geometry 
eigenstates. The implementation of this is currently a work in progress by the authors. 
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Lastly, the duality between general relativity and shape dynamics is not accidental. There is 
a general mechanism at work that can be seen in a variety of classical gauge theories. This is a 
consequence of the second result of this paper: given a gauge fixing that satisfies certain inegrability 
conditions, one can construct a dual theory by replacing the original gauge generators with the gauge 
fixing conditions. 
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A Direct phase space geometry of exchange 

We now give a direct and formal proof of the steps involved in the dualization in geometro dynamics. 
To justify the use of geometrical arguments, we should have that the regular value set 

nite)- 1 ^)}, 

i 

for our constraints Xii should initially form a first class submanifold of phase space. In the infinite- 
dimensional context, the immersions need to be splitting. This is a non-trivial condition that does 
not follow solely from the finite-dimensional result of manifolds formed from regular points. For 
definiteness, let us take the constraints to be T~L and H a . For these, one can (almost 18 ) show that 
the intersection of the regular points is indeed a manifold using the Fredholm alternative [21, 22]. 
Now we are ready to state our process in a geometrical setting. 

A.l Formal Steps 

1. Extend ambient phase space 

We define the extended ambient manifold, Te := T x T*C, by using a global parametrization 
of the additional coordinates through ((f), tt^). Restricting to the subspace parametrized by (j) 
(which is nonetheless more conveniently parametrized, albeit redundantly, by <j>) we find that 
(n^) has trivial action on all of the variables. Thus, we can formally restrict our attention to 
the space 

r e := r x T*(C/V) (92) 

conveniently (and redundantly) parametrized by (gij, ir l i , eft, ir^). Note that the redundancy 
has a nontrivial dependence on T. 

2. Stiickelberg Mechanism 

18 For the proof that this intersection is indeed a manifold, one needs to further use the constraint that tt = c(i), i.e. 
that the mean extrinsic curvature be a spatial constant. This is precisely the condition that arises in our procedure. 
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For T : T — > C°°(S), we define the natural extension 

TT : T e -> C°°(E), (93) 

This is the simple transformation TJ r (gij,ir i:! ) = ^(T^j, T^tt*- 7 ). Under this extension, we 
have the following identities: 

{V(8),Tg ab } = {ir lj> (6),Tg ab } and {V(6), Tir ab } = {^(9), Tvr afe }, (94) 

where 9{x) is a smearing function. Thus, we have the following quantity acting trivially on 
the image of T: 

C(x) = tt(x) - y/g(x)(ir) g - n^x). (95) 

Imposing this constraint on functions on r e is an efficient way to require our dynamical 
systems to be in a one to one relation with the one defined in the original T. This constraint, 
arising from the Stiickelberg mechanism, eliminates the extra degrees of freedom that have 
been introduced. Because of our redundant parametrization, C(x) contains the redundant 
constraint (C) = 0. 

3. Impose tt^ = and separate first and second class parts. 

it j, = trivially forms a first class submanifold of r e . Its symplectic flow is equivalent to that 
of T> on the image of T and it is thus tangent the two first class submanifolds C(x) ~ and 
T^U a « 0: 

{C,vr4 = and {TH a {N a )^}^Q. (96) 

Their intersection consists of the lift through the T map of the space of transverse momenta 
with spatially constant trace given by (it). This can be shown to form a manifold. 

Regarding the space generated by the remaining constraints T%(x) (one at each point), there 
is still one direction, among the infinitely many, that is tangent to 

nWOzO}- 1 ^) f]{TH{y)}-\<S). (97) 

This is 

{TS(N ),^(x)} oc TS(x) w (98) 

where No[gij; (f>, x). Thus, the direction given by H g f = J (fix Nq(x)TS(x) is first class. 
Furthermore, (98) has a homogeneous term that includes the integral of N so that its solution 
is only determined up to a constant scaling: i.e. if iVo is a solution so is ciNq for a G R/{0}. 
We, therefore, have existence and uniqueness of solutions only on the domain C°°(S)/V. We 
can split TS(x) into a part that is first class with respect to all the constraints and a part 
that is second class: 

TS(N ) and TS(x) := TS(x) - TS(N ), (99) 

respectively. 
• Intemezzo 

We now show how these properties arise. The homogeneous lapse fixing equation resulting 
from the Poisson bracket {H(N), 7r(x)} = is 

(V 2 + f(x))N(x) := AN(x) = 0, (100) 
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where we define the second order elliptic self-adjoint differential operator A = (V 2 +/(x)) and 
f(x) = —(1 + 22) in the case we are considering. There are many different ways to prove that, 
as a differential operator, A is an isomorphism from smooth functions to smooth functions 
(see for example [19] or [23] for the construction of the respective Green's functions 19 ). Taking 
this to be the case, A is an invertible operator to which one can find solutions for appropriate 
boundary data. 

The inhomogeneous analogue to (100) is 

AN(x) - (A(N)) = {H(N), (vr - (ir)Jg)(x)} = 0. (101) 

Using the above arguments, we have a one parameter family of solutions to (101), that is, 
whenever AN(x) — (A(N)) = c where c is a given non-zero constant. Suppose we have a 
solution, No(x,gij], such that c = 1. Applying the canonical transformation T, we get 

(TA)N(x) - ((TA)N) = {TH(N), (tt - (tt)^)(x)} = T{H(N), (vr - (tt)^)(x)}. (102) 

A solution, Nq(x; gij,(p], to (102) is 

N{xig ij) <j>]=T^N(x,g ij ], (103) 

where N(x,gij] is the solution to the untransformed equation (101). 

• End of intermezzo 

4. Solving the constraints 

For the unconstrained conformal transformations, which we will denote by T, 

TS{x) :T xT*{C) ^C, (104) 

where we used that C°°(E) ~ C. Since these equations do not depend on tt^, we can fix 
^4>{x) = f(x). Then 

TS(x)^ =f{x) : T x C7°°(S) -> C°°(E). (105) 

Denoting the derivative in the second coordinate, the one parametrized by <j>, by a subscript 
C, we have, from equation (100), that the linear operator 

5 C TS := = {TH(x),^(y)} = A(x)6(x - y) (106) 

is a topological linear isomorphism 20 between the spaces C°°(£) and C°°(S). However, for 
the volume preserving conformal transformations 

TS{x)^ =m : r x C7°°(5])/V -> C°°(S), (107) 

we explicitly derived, from properties of the full conformal transformations, that the linear 
operator 

ScTS := = {TH(x),^(y)} (108) 

does not yield a topological linear isomorphism between the spaces C°°(£)/V and C°°(S). 



19 This can also be extended to distributional domains. 



20 We assume that, for all practical purposes, we can carry on as if these were Banach spaces. 
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The adjoint "matrix" has a one dimensional kernel given by TH(Nq). However 

TS{x)^ =f(x) : T x C°°(S)/V -> C°°(£)/V 

has the infinitesimal transformation 

5TV(x) - — 

ScTS := -^±± = {TH{x)^(y)}. (109) 

Having taken out the kernel, this is necessarily an invertible operator. Not only can we form 
the Dirac bracket using {TH(x), ^(y)} -1 but we can now use the implicit function theorem 
for Banach spaces to assert that 

Proposition 2 There exists a unique 4>q : T — > C/V such that 

(fsrHo) = {{g lv TT l \Ua^y 3 l^) I (9v^ ij ) e r}. 

5. Construct the theory on the constraint surface 

We have then a surface in T e , defined by ir^ = and <fi = fa, on which T% = 0, and whose 
intrinsic coordinates are gij,TT lJ . Furthermore, the Dirac bracket on the surface exists and we 
can use it to help prove theorems on the constraint surface, which now has the symplectic 
structure 

{') '}\reduced -~ {'; '}db \<p=<f>o, n$=0 ~ {' \4>=<f>o ,n $=0 1 ' |0=</>o,ti>=o} • (HO) 

One can immediately see from (110) that the first class constraints D,TT-L a and (TNoH) 
remain first class. Furthermore, we can verify that 

{•,T(fX}\reduced = Tffoi^X} (HI) 

for the remaining first class constraints \ ( as they are represented in the original phase space). 
Thus, T^ Q is a canonical transformation for any of the first class constraints, giving us all we 
need for the dynamics. Finally, we reconfirm the conclusions of the main text: that we have 
a dual Hamiltonian 
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#duaJ = N{N U) + I d 3 x (A(s) (vr(x) - (vr) y/g) + p a {x)H a {x)) 
Jt, 

in r with the first class constraints 

(N H), V,U a (113) 
and the physical configuration space A^/(Diff(S) x (C/V)). 

• Final remark 

We could have used the operator 

(TA)S(x, y) = {TH(x), 7r(y)} = {TH(x), ^(y)}, (114) 

which generates the full group of conformal transformations, instead of T. This operator 
is invertible because of the invertibility of A. Also, we could have both used the implicit 
function theorem and formed a Dirac bracket. However, the homogeneous equation only 
has the solution N = for compact manifolds without boundary. These don't generate any 
dynamics. 
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